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Abstract-  High  intracellular  calcium  conditions  cause  a 
calcium-activated  transient  inward  current  ( lu )  that  can 
provoke  oscillations  in  membrane  potential  called  delayed 
afterdepolarizations  ( DAD ).  The  current  I,i  comprises  of  the 
sodium-calcium  exchange  current  ( Ipjaca )  and  the  calcium- 
activated  chloride  current  ( Ici-ca  )•  Lindhiad,  Murphey,  Clark 
and  Giles  developed  a  mathematical  model  (  LMCG  model )  of 
the  rahhit  atrial  AP.  In  this  study,  a  modified  AP  LMCG  model 
that  includes  Ici-ca  is  used  to  evaluate  the  contrihution  of  Ici-ca  to 
develope  DADs.  Our  results  suggest  that  although  iN^ca  is  the 
main  component  of  lu  (  65  %  ),  Ici-ca  niay  play  a  significant  role 
in  DAD  generation.  Even  more,  the  Ici-ca  blockade  could  inhibit 
the  DAD  propagation  and  trigger  activity  associated  to  high 
[Ca^^i  condition,  in  atrial  tissue. 
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1.  Introduction 

Delayed  afterdepolarizations  (  DADs  )  are  oscillations  in 
membrane  potential  ocurring  after  completion  of  the  action 
potential  (AP).  DADs  are  an  important  mechanism  for 

cardiac  arrhythmias  [1,  2].  They  usually  occur  under 
conditions  in  which  [Ca^"^]!  appears  to  be  high.  The 
mechanism  underlying  DAD  is  a  transient  inward  current 
( Iti )  activated  by  spontaneous  Ca^"^  release  from 

sarcoplasmic  reticulum. 

Although  the  ionic  nature  of  1,1  is  still  subject  to  debate, 
recent  experiments  in  ventricular  and  purkinje  myocytes 
indicate  that  I,;  is  composed  by  an  electrogenic 

sodium-calcium  exchange  I^^ca  and  a  calcium-activated 
chloride  current  Ici-ca  [!]•  These  studies  suggest  that  the  main 
component  of  !„  is  Ij,jaCa  (  about  60-^80  %  )  but  the  role  of 
Ici-ca  could  be  important. 

Ici-ca  is  also  known  as  the  second  component  (  1,02  )  of  a 
transient  outward  current  (  I,o  )  that  is  present  in 

depolarization.  Ici-ca  magnitude  is  lower  than  Itoi,  a  potasium 
current,  at  normal  conditions  and  slow  rates.  It  is  dependent 
on  [Ca^"^]!  ,  time-independent  and  it  presents  a  voltage- 
dependent  outward  rectification  [4,  5,  6,  7].  1,02  is  responsible 
for  the  maintenance  of  AP  duration  at  high  rates  [5]. 

We  developped  a  comprehensive  model  for  the  Ici-ca ,  that 
was  presented  in  a  previous  comunication  [8]. 

This  model  has  been  introduced  in  the  Lindblad,  Murphey, 
Clark  and  Giles  mathematical  model  (  LMCG  model )  [3]  of 
the  rabbit  atrial  myocyte  that  reproduces  the  AP  in  a  variety 
of  conditions. 

The  present  paper  pretends  to  clarify  the  role  of  Ici-ca  in 
DAD  generation  and  in  the  induction  of  triggered  activity. 


II.  Methodology 


A  set  of  mathematical  equations  describes  the  dependence 
of  Ici.ca  on  intracellular  Ca^"^,  as  well  as  on  ionic  concentrations. 
The  current  through  Cl-Ca  channels  is  described  as  a  product  of 
three  terms:  first  term  is  expresed  using  the  Goldmann- 
Hodgkin-Katz  (GHK)  equation,  the  second  term  ( fca )  is  the 
dependence  on  calcium  and  follows  a  Hill-type  equation,  and 
the  last  term  (  /?c  )  is  the  rectification  exhibited  by  Cl-Ca 
channels,  which  is  due  to  a  voltage-dependent  blockade  caused 
by  intracellular  cations  (  Table  1 ) . 

Table i 

CALCIUM-DEPENDENT  CHLORIDE  CURRENT  EQUATIONS 
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Model  Parameters: 

Kmcrca  =  150,210'^  mM/1  pci  =  1,1712-10'\S 


The  mathematical  description  of  Ici.ca  was  included  into  the 
atrial  AP  LMCG  model.  Software  programs  were  written  in 
ACSL  language  using  Gear  stiff  algorithm  to  solve  the 
nonlinear  system  of  differential  equations  that  results  from 
the  AP  myocyte  model.  For  studying  propagation,  the  tissue 
is  modeled  as  a  one-dimensional  1-D  segment  (  500  cells  ). 


V(mV) 


Ligure  1.  Effect  of  varying  Ici-ca  conductance  (10,  25,  50,  75,  100,  150, 
200,  250,  300  %  of  nominal )  on  model  action  potential  morphology 
during  stimulation  at  2  Hz. 
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The  1-D  model  was  programmed  in  C++  language,  the 
cable  equation  was  solved  using  a  central  difference  scheme 
in  space  and  Crank-Nicholson  method  in  time.  Both,  isolated 
myocyte  and  multicellular  models,  were  implemented  in  a  PC 
machine. 

Protocol  stimulation  follows  the  nominal  values  of  LMCG 
model  for  stimulation  current  (Is,)  at  2  Hz  rate  [3]. 


III.  Results 

Figure  1  shows  the  effect  of  varying  the  magnitude  of  the 
nominal  Ici.ca  conductance  on  action  potential.  Blockade  of 
Ici.ca  increase  the  action  potential  duration. 

Figure  2  shows  the  time  course  of  different  membrane 
currents  during  the  AP.  Under  normal  conditions,  Ici-ca  is 
present  mainly  in  depolarization. 

Figure  3  shows  the  parameters  fca  and  Rc  used  to  model 
Icica-  In  order  to  compare  the  model  Ici-ca  I-V  plots  with 
experimental  data  [5,  6],  a  train  of  voltage  steps  of  100  ms 
from  -  70  to  V  was  simulated  at  0.1,  1  and  2  Hz.  Plot  C 
represents  the  normalized  model  results  of  peak  current  vs. 
experimental  normalized  data. 

In  order  to  study  !„  we  considered  DAD-generating 
conditions  as  intracellular  calcium-overload.  Calcium- 
overload  is  induced  by  an  initial  [Ca^'^Jo  higher  than  nominal. 

In  figures  4,  5  and  6,  a  train  of  10  AP  was  stimulated  but 
only  the  last  AP  was  represented. 

Figure  4  shows  the  membrane  potencial  at  different  initial 
[Ca^'^lo  (  A;  4,  B:  5,  C:  6  and  D;7  mM/1)  in  control  simulation 
(a)  and  when  Ici-ca  (b)^  iNaCa  (c)  and  both  (d)  are  forced  to 
zero  at  4.8  s. 

It  can  be  observed  that  at  higher  calcium-overload 
conditions  the  delay  between  AP  and  DAD  is  lower  (t4<  t,). 

For  [Ca^'^lo  of  4  and  5  mM/1,  the  inhibition  of  Ici-ca  has  not 
significant  effect  on  DAD  amplitude  (  differences  are  below 
than  11%).  However,  for  [Ca^'^Jo  of  6  mM/1  the  inhibition  of 
Ici-ca  provokes  a  reduction  of  58  %  in  DAD  amplitude. 

Figure  5  shows  the  case  of  an  initial  [Ca^'^]o=  6  mM/1.  In 
plot  A  can  be  observed  the  effect  of  Ici-ca  blockade  at  5.2  s 


Figure  3.  Parameters  used  to  model  calcium-activated  chloride  current  ( Ici-ca  )•  A:  calcium-dependent  open  channels  fraction.  B:  voltage-dependent 

blockade.  C:  I-V  normalized  relation.  Data  points  represent  experimental  data. 


Figure  4.  A  secuence  of  10  PA  at  2  FIz  was  made.  Calcium-overload  is  provoked  by  increasing  [Ca*^]o  to  4,  5,  6  and  7  mM/1  (  nominal  value  is  2.5  mM/1 ).  At 
4.8  s  we  force  to  zero  Icica,  iNaCa  and  both  ( traces  b,  c,  d  )  comparing  them  with  control  ( trace  a ).  The  results  in  plot  A  and  B  show  similar  magnitude  in 
traces  a  and  b.  The  plot  D  shows  a  great  attenuation  of  DAD  traces.  Only  plot  C  seems  to  have  a  relevant  difference  between  DAD  amplitudes  (traces  a,  b  ). 


(b)  vs.  control  (a)  on  DAD.  In  plot  B  can  be  observed  that  the 
contribution  on  of  Ici-ca  ^NaCa  is  35  %  and  65  % 
respectivelly. 

To  study  DAD  propagation  we  use  the  1-D  model  in  a 
segment  of  500  cells. 

Figure  6  shows  the  membrane  potential  in  the  cell  number 
1,  125,  250,  375  and  500.  As  it  can  be  observed,  in  plot  A  the 
DAD  that  appears  in  cell  #1  is  propagated  like  an  AP  in  the  cell 
#500.  However,  when  laca  is  blocked,  plot  B,  the  oscillation  in 
the  membrane  potential  of  cell  #1  is  not  propagated  like  an  AP 
in  cell  #500. 

IV.  Discussion 

The  main  goal  of  the  present  study  is  to  find  out  the  role  of 
Icica  as  a  component  of  on  DAD  formation  and  propagation 
as  AP  in  the  atrial  tissue. 


mV 


A  mathematical  model  of  the  Icica  has  been  developed. 
Whenever  possible,  we  use  data  directly  measured  on  rabbit 
atrial  cells  to  derive  model  parameters  [5,  7].  When  rabbit  atrial 
data  were  not  available  to  completely  caracterize  the  current  we 
use  reliable  data  of  dog  and  rabbit  ventricular  cells  [4,  6]. 

The  ionic  nature  of  is  still  subject  to  debate.  Several 
authors  propose  two  ionic  currents  to  contribute  to  1^; 
electrogenic  I^aca  and  Igca  ,  but  it  is  not  clear  in  which 
percentage. 

Recent  data  from  Zygmunt  el  al.  [1]  indicate  that  I^aca 
component  is  the  main  responsible,  representing  about  60%  of 
the  total  calcium-activated  current  at  the  resting  potentials  ( Iti ). 
Other  authors  like  Verker  et  al.  [2]  give  values  near  to  80%. 

Our  results  show  that  Ici-ca  takes  values  near  to  35% 
and  iNa-ca  takes  values  near  to  65%  of  the  total  calcium- 
activated  current  at  the  resting  potentials  ( lu  ). 


Figure  5.  DAD  after  10  PA  @  2  Flz,  [Ca^*]o  =  6mM/l,  CTRL  and  lcica=0.  A:  DAD  plot.  B:  la  CTRL  composition  %. 
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Figure  6.  A  secuence  of  10  PA  at  2  Hz  in  a  1-D  segment  of  500  cells  was 
A  show  the  tissue  propagation  of  the  DAD.  In  the  plot  B  we  force  to  zero 
lower  than  the  subthreshold  level  for  triggered  AP. 


mV 


de.  Calcium-overload  is  provoked  increasing  [Ca*^]o  to  7  mM/1.  The  plot 
:a  at  t>4.6  s  and  AP  <  10  %  of  maximum  amplitude;  in  this  case,  DAD  is 


V.  Conclusion 

Our  results  suggest  that  Ici-ca  could  play  an  important  role  in 
DAD  generation  mechanisms.  The  blockade  of  Ici-ca  can  reduce 
DAD  amplitude  to  prevent  DADs  to  reach  the  trigger  threshold 
of  action  potentials. 
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